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While strong topological insulators (STI) have been experimentally realized soon after their the¬ 
oretical predictions, a weak topological insulator (WTI) has yet to be unambiguously confirmed. A 
major obstacle is the lack of distinct natural cleavage surfaces to test the surface selective hallmark 
of WTI. With a new scheme, we discover that Bi 4 X 4 (X=Br, I), stable or synthesized before, can be 
WTI with two natural cleavage surfaces, where two anisotropic Dirac cones stabilize and annihilate, 
respectively. We further find four surface state Lifshitz transitions under charge doping and two bulk 
topological phase transitions under uniaxial strain. Near the WTI-STI transition, there emerges a 
novel Weyl semimetal phase, in which the Fermi arcs generically appear at both cleavage surfaces 
whereas the Fermi circle only appears at one selected surface. 


The discovery of topological insulators [1-3] has led to 
an ongoing revolution deepening our fundamental under¬ 
standing of quantum materials. The controllable bulk 
topological quantum phase transitions and the protected 
spin-momentum locked surface states may ultimately 
lead to unprecedented advances in technologies, e.g., the 
Majorana-based fault-tolerant quantum computing [4]. 
Historically, strong topological insulators (STI) and weak 
topological insulators (WTI) were predicted together to 
exist in three dimensions [5-7]. Astonishingly, STI are 
not at all rare in nature. Over a dozen documented ma¬ 
terials have been identified as STI under readily acces¬ 
sible experimental conditions [8-12]. Notably, their lat¬ 
eral structure offers a superior advantage to observe the 
hallmark of STI, a surface Dirac cone, without intended 
surface passivation that could be challenging. However, 
the unambiguous experimental confirmation of WTI is 
still elusive [12]. Thus far there have been two routes 
to construct WTI, i.e., to stack weakly coupled layers of 
quantum spin Hall insulators (QSHI) [13, 14], or to en¬ 
gineer a superlattice of alternating layers with multiple 
band inversions [15, 16]. Evidently, these designer WTI 
pose extreme experimental challenges. In the first route, 
the protected surface metallicity will only be present at 
non-cleavage surfaces (non-parallel to the layers), yet sur¬ 
face roughness and dangling bonds would prevent us from 
observing this hallmark [14]. In the second route, though 
the surface states can survive at the cleavage surface (par¬ 
allel to the layers), the intra- and inter-layer couplings 
within a supercell must be fine tuned in material synthe¬ 
sis. Therefore, to unambiguously determine the existence 
of WTI and to explore the exotic phenomena uniquely 
hosted by WTI, it is crucial to develop a new route. 

We propose to realize the WTI in a van der Waals 
(vdW) material that is a periodic stack of one dimen¬ 
sional atomic chains. Such a WTI possesses two natu¬ 
ral cleavage planes, enabling respective observations of 
stabilization and annihilation of distinct surface Dirac 
cones. We discover that /3-Bi4X4 (X=I, Br) can be such 
a class of WTI. Markedly, with a chemical formula as 
simple as the prototype STI Bi 2 X 3 (X=Se, Te), /3-Bi4X4 


(X=I, Br) are real crystalline solids rather than artihcial 
periodic heterostructures. /3-Bi4l4 has been successfully 
grown as a large crystal before [17-20], and /3-Bi4Br4 is 
demonstrated to be similarly stable [21]. After revealing 
their unique crystal structure, topological band proper¬ 
ties, and surface Lifshitz transitions, we will further ex¬ 
amine their topological phase transitions under uniax¬ 
ial strain. Intriguingly, a novel Weyl semimetal (WSM) 
phase emerges near the WTI-STI transition; in addition 
to the existence of Fermi arcs at both cleavage surfaces 
as the case of WSM, a Fermi circle only exists at one 
selected surface reminiscent of the case of WTI. As pro¬ 
totype WTI, ,0-Bi4X4 (X=I, Br) offers a new platform 
for exploring exotic physics with simple chemistry. 



FIG. 1 : Crystal structure and Brillouin zone of / 3 - 
314X4 (X=I, Br). (a) The crystal structure of the con¬ 
ventional cell of P-BiPAi. Biin (Biex) and Bi[n (Bi(,x) atoms 
are related by inversion symmetry, (b) The first Brillouin 
Zone of /3-Bi4X4 with eight time-reversal-invariant momenta. 
G are the reciprocal lattice vectors. 
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FIG. 2: Band structures and surface states of /3-Bi4X4 (X=I, Br). (a)-(b) The bulk band structures for jS-Huli without 
and with the SOC. The size of red (blue) dots indicates the weight of the relevant orbital of Biin (Biex) atoms; the symbols 
± label the parities of the bands at L and M points; the dashed lines are the Fermi levels, (c) The (001) cleave-surface states 
of /S-BuB. (d)-(e) The same as (a)-(b) but for /?-Bi 4 Br 4 . (f) The (100) cleave-surface states of /?-Bi 4 Br 4 and (g) their Lifshitz 
transitions. Vertical dimensions in (g) are magnified by 33 times for clarity, (h) The illustration of the stacking orders in a- 
and /3-Bi4Br4. The P (a) phase is a stack of single (double) (001) layers; each layer denoted by a green line is a QSHI. 


Topological Band Structures 

Both Bi 4 l 4 and Bi 4 Br 4 have stable a and /3 phases. 
The two phases crystallize in the same monoclinic space 
group Cf/j {C2lm), and differ only in the way their build¬ 
ing blocks are stacked [17-20]. The a phases turn out to 
be trivial insulators [21]. As the /3 phases have rather 
similar lattice structures, by way of illustration we re¬ 
fer to /3-Bi4l4 that has been experimentally synthesized, 
and the explicit lattice constants and atom positions of 
/3-Bi4Br4 are given in Table SI [21]. Figure la shows 
the conventional cell of jP-Biili, in which a = 14.386 A, 
b = 4.430 A, c = 10.493 A, and /3 = 107.87°. One unit 
cell consists of four I atoms and four Bi atoms that can 
be divided into two types. The two internal Bi atoms, 
labeled as Biin and BiJ^^, form zigzag atomic chains with 
the nearest neighbor distance of 3.04 A. The two exter¬ 
nal Bi atoms, labeled as Biex and Bi^^, are each bonded 
to four I atoms with the Bi-I distance of 3.14 A and 
to one internal Bi atom with a distance of 3.06 A. The 
crystal has two independent symmetries, spatial inver¬ 
sion and mirror reflection. Under (010) mirror reflection, 
the four Bi atoms are invariant. Upon spatial inversion, 
the two Bi atoms are interchanged for each type. As sug¬ 
gested by Fig. la, /S-BUU is a periodic stack of atomic 
chains aligned to the b direction. Indeed, we find that 
the interlayer binding energies for (100) and (001) planes 
are about 20 meV/A^, comparable to 12 meV/A^ for 
graphite [22] and 26 meV/A^ for M 0 S 2 [23]. Evidently, 
/3-Bi4X4 (X=I, Br) are typical van der Waals (vdW) ma¬ 


terials, but with two cleavage surfaces. 

We employ the HSE hybrid functional method [24], 
which is believed to be more sophisticated and accu¬ 
rate than the GGA and EDA methods, to carry out 
our density functional theory (DFT) calculations of the 
electronic band structures of /S-BUU and /3-Bi4Br4. The 
vdW corrections [24] and the lattice relaxations are also 
taken into account to optimize the crystal structures. 
The projected band structures of /S-BUU without and 
with SOC are shown in Figs. 2a-2b, respectively. The 
bands near the gap are mainly contributed from the Bi 
Px orbitals. Without the SOC, the conduction and va¬ 
lence band edges at both L and M points are respectively 
from the Biex and Biin orbitals, which are identified to 
exhibit opposite parities. When the SOC is included, the 
constituents and parities of the conduction and valence 
bands remain the same at the L point, whereas they are 
inverted at the M point. Based on the Fu-Kane parity 
criterion [25], we can conclude that /S-BUU is a STI, with 
(1; 110) Z 2 invariants and a 39 meV indirect gap. As 
plotted in Fig. 2c, we further obtain one surface Dirac 
cone at the (001) cleavage surface, which is consistent 
with the conclusion that /S-BUU is a STI. Notably, the 
Dirac cone is highly anisotropic, because the crystal only 
has inversion and reflection symmetries, and because the 
surface is parallel to the atomic chains. 

We now study the band structure of /3-Bi4Br4, which 
share evident similarities to the case of /3-Bi4l4. How¬ 
ever, the band inversion occurs at both M and L points 
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in / 3 -Bi 4 Br 4 , as shown in Figs. 2d-2e. It follows from 
the Fu-Kane parity criterion [25] that /I-Bi 4 Br 4 is a WTI 
with (0; 001) Z 2 invariants and a 32 meV indirect gap. 
We further calculate the surface states for the two natural 
cleave surfaces, i.e., the (100) and (001) surfaces. Since 
M and L are projected into the same point at the (001) 
surface Brillouin zone, the (001) surface are anticipated 
to host neither protected surface states nor dangling bond 
states. Our calculations verify this physical picture (not 
shown). In contrast, at the (100) surface, M and L are 
projected into two distinct points G and A. As antici¬ 
pated and shown in Fig. 2f, our calculations identify two 
(connected) surface Dirac cones at G and A points; the 
dangling bond states are completely absent. 

The (100) surface states of / 3 -Bi 4 Br 4 also exhibit 
prominent anisotropy, for the reasons that we have dis¬ 
cussed above. The direction with the much larger Fermi 
velocity is along 6 , in which the atomic chains are aligned. 
Moreover, the two surface Dirac points at G and A are 
at different energies. These two features enrich the Fermi 
surface topology of the (100) surface states. Fig. 2g char¬ 
acterizes the four corresponding Lifshitz transitions, as 
we now explain in ascending energy. At zero energy, the 
two hole pockets around G and A are connected. As the 
energy increases, e.g., to 3 meV, the two pockets are dis¬ 
connected after the first Lifshitz transition. At 7 meV, 
the second transition occurs; the hole pocket at A con¬ 
tracts into a Dirac point, followed by the emergence of 
an electron pocket. At higher energies, the disconnected 
electron and hole pockets coexist up to the third tran¬ 
sition at about 24 meV, in which the hole pocket at G 
contracts into a Dirac point and then an electron pocket 
emerges. Ultimately, the fourth transition occurs and the 
two electron pockets become connected. 

We note by passing that, in sharp contrast to the 
case of /3 phase, both a-BUU and a-Bi 4 Br 4 are normal 
insulators (NI) [21]. The contrasting topological prop¬ 
erties of the a and (3 phases can be understood in the 
way illustrated by Fig. 2h. The (3 phase is a stack of 
single ( 001 ) layers, whereas the a phase is a stack of 
double (001) layers. For Bi 4 Br 4 , the single (001) layer 
was demonstrated to be a QSHI with a 0.18 eV gap [23]. 
Given the weak interlayer couplings, the (3 phase is a 
WTI with similar band inversions at = 0 and = tt 
planes, whereas the a phase is trivial since each bilayer 
is a Nl. For BUU, the single (001) layer was shown to 
be close to the QSHl-Nl critical point [23]. Thus, the a 
phase is definitely a Nl, yet the /3 phase is not necessarily 
a WTI, depending on the details. However, as we will 
show, / 3 -Bi 4 l 4 can become WTI under uniaxial strain. 

Topological Phase Transitions 

We now construct an effective model for the /3 phases. 
The time-reversal-invariant M and L points have the lit¬ 
tle group G|^ with two independent symmetries, spatial 
inversion V = Tz and mirror reflection A4y = iay. Here 
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FIG. 3: Topological phase transitions and band inver¬ 
sions under uniaxial strain along the a axis. The phase 
diagrams versus the uniaxial strain are depicted in (a) for /?- 
BUU and (b) for d-Bi 4 Br 4 . The corresponding direct gaps 
at L and M points are shown in (c) for / 3 -Bi 4 l 4 and (d) for 
/l-Bi 4 Br 4 . A solid (open) dot denotes a direct (indirect) gap; 
the minus sign of a gap indicates a band inversion. 


T and (T are the orbital and spin Pauli matrices. By con¬ 
vention we choose the time-reversal operator T = iKay 
with K the complex conjugation. Given the three sym¬ 
metries, to the linear order the k ■ p Hamiltonians near L 
and M may be written as [29] 

T-C = vlk^ayTx + vlkyOxTx -f vlk^Ty + mGz + c\ (1) 

where i refers to the M or L point, v are the velocities, 
c is the energy offset, and m is the energy gap with m < 
0 denoting an inverted gap. For j3-QiiG only is 
negative; for /?-Bi 4 Br 4 both and are negative. 

With Eq. (1), we elaborate the strain effects on (3- 
Bi 4 X 4 (X=I, Br). The strain tensor is a rank-2 tensor 
and invariant under time reversal. With these two facts 
and the original symmetries of %, the strain induced per¬ 
turbations to the lowest order take the form of 


5T-C — ST-Lq + (eiiA 44 -I- 622^22 + £^ 33^33 -I- ei3A43) Tz, (2) 

where A^- are deformation potentials, and SHq is a rigid 
shift of all bands with a form similar to that in the paren¬ 
theses. Clearly, applying strain along the a, 6, or c axis 
would change the direct gaps and hence the band inver¬ 
sions at the M and L points. Moreover, the direct gaps 
should exhibit linear dependence on the small strain, and 
the inversion symmetry is respected in the deformation. 
It follows that the topological phase transitions among 
NI, STI, and WTI can be tuned by the uniaxial strain in 
/ 3 -Bi 4 X 4 (X=I, Br). 

This prediction can be verified by our DFT calcu¬ 
lations [24] with strain along the a axis, as shown in 
Fig. 3. In this case, the deformation potentials at M and 
L points exhibit the same sign and similar magnitudes. 
For the STI /S-BUU, under more than 1.5% tensile 
strain, the inverted bands at M become un-inverted 
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FIG. 4: Phase diagram and surface states for compos¬ 
ite Weyl semimetal, (a) Schematic phase diagram of WTI, 
CWSM, STI, WSM, and NI. The two curves depict the direct 
gaps near the relevant L and M points, (b) Existence of two 
open Fermi arcs at the (001) surface of the calculated CWSM. 
(c) Coexistence of two open Fermi arcs and one closed Fermi 
circle at the (100) surface of the calculated CWSM. In (b) 
and (c), the dimensions have been magnified for clarity. 

while those at L remain not inverted, yielding a STI to 
NI transition. On the other hand, under more than 3.3% 
compressive strain, the bands at L become inverted 
while those at M remain inverted, producing a STI to 
WTI transition. For the WTI /I-Bi 4 Br 4 , the inverted 
gaps at L and M both increase with increasing the 
compressive strain, resulting a larger gap WTI. As the 
tensile strain increases, however, there are two successive 
transitions from the WTI to a STI and then to a NI, 
because the bands at L become un-inverted prior to 
those at M. Figures 3a-3d plot the phase diagrams 
and the direct gaps at L and M, which exhibit linear 
dependence on the strain predicted by our effective 
model. Similar phase transitions can also be driven by 
strain along other axes [21]. These results are suggestive 
of a feasible way to engineer or stabilize the WTI phase. 

Composite Weyl Semimetals 

In the presence of both V and T symmetries, all bands 
are doubly degenerate and a Dirac point accounts for the 
NI-STI trainsition. The Dirac point is unstable unless 
there is a symmetry protection [30-33]. However, when 
V or 'T symmetry is broken, the Dirac point splits into 
pairs of Weyl points, each of which is locally protected 
by the Chern number of a constant-energy surface 
enclosing it. One might naively think the same WSM 
phase [34-39] naturally emerges near the WTI-STI 
transition. In fact, the emergent phase turns out to be 
novel. We now examine such a phase using /I-Bi 4 l 4 under 
3.3% compressive strain (see Fig. 3a). The V symmetry 


can be further broken by a possible energy difference 
between the Bi and Bi' atoms. Because of the lack of 
a rotational symmetry, only two pairs of Weyl points 
appear at the WTI-STI transition; the bulk is a standard 
WSM. At the (001) surface, there are two open Fermi 
arcs connecting the four projected Weyl nodes, as shown 
in Fig. 4b. At the (100) surface, surprisingly, there exists 
one closed Fermi circle in addition to the anticipated 
Fermi arcs, as shown in Fig. 4c. We dub this novel phase 
the composite WSM (CWSM). When V asymmetry is 
small, while the Weyl points emerge near L, the band 
gap remains inverted at M. As a consequence, there 
are surface states locally, separated by a large crystal 
momentum, at the (100) surface where L and M project 
into distinct points A and G. In contrast, at the (001) 
surface where L and M project into the same point G, 
the Fermi arcs appear attribute to the nontrivial Chern 
numbers, whereas the Fermi circle disappears because 
of its strong scattering with the metallic background. 
Intriguingly for a CWSM, the Fermi arcs generically 
appear at any surface like the case for WSM, whereas 
the Fermi circle only appears in selected surfaces like 
the case for WTI. In the most generic topological phase 
diagram, as sketched in Fig. 4a, the novel CWSM phase 
emerges near the WTI-STI transition when V or T 
symmetry is broken. 

Discussions 

One may wonder whether ,5-Bi4l4 and /3-Bi4Br4 are 
TCIs [26], since a (010) mirror symmetry is present. In 
a mirror invariant plane, the two spin subspaces, related 
by T symmetry, decouple and have opposite Chern num¬ 
bers ±A^. There must be N pairs of protected Fermi 
points in the mirror invariant line (MIL) at the surface, 
independent of the Fermi level. Similar to Bi 2 Se 3 [27], 
the STI /3-Bi4l4 is also a TCI {N = 1), because there are 
always two surface-state Fermi points in the MIL GX, 
as shown in Fig. 2(c). In contrast, the WTI /3-Bi4Br4 is 
not a TCI {N = 0), because Fermi points are not nec¬ 
essarily present along the MIL GA, as seen in Figs. 2(f) 
and 2(g). These arguments are consistent with our direct 
calculations on the mirror Chern numbers [28]. 

,d-Bi 4 X 4 (X=I, Br) can be prototype WTI and offer 
a new platform for exploring exotic physics with simple 
chemistry. In addition to the strain-induced topological 
phase transitions and CWSM phase, the two (100) sur¬ 
face states is particularly appealing, e.g., the spin tex¬ 
ture, the Landau level crossing, and the exciton conden¬ 
sation triggered by the nesting of the electron and hole 
pockets. In the thin-film QSHI limit, the small velocity 
along the a or c axis can be utilized to study the helical 
Luttinger liquid [40] and the Z 4 parafermions [41-43]. As 
for the CWSM, physics can be enriched by single-particle 
couplings or many-body interactions between the open 
Fermi arcs and the closed Fermi circle. 

Finally, WTI have been proved to be strong against 
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disorder [44-46], as long as U{V), T, and translational 
symmetries are respected on average. Breaking any 
symmetry at the (100) surface of /3-Bi4Br4 or strained (3- 
Bi 4 l 4 may induce an exotic phenomenon. A topological 
defect, e.g., a screw dislocation [47] or a step edge [48], 
can break a translational symmetry binding a helical 
edge modes. A Zeeman field can break T symmetry 
yielding a quantum anomalous Hall effect with a Chern 
number between —2 to 2, which is tunable by the field 
orientation [27]. A proximity coupling to an sj- wave 
(e.g., iron-based) superconductor, can break U{1) gauge 
symmetry producing a Z 2 topological superconductor 
with a Majorana Kramers pair [49]. 

Computational Methods 

Our DFT calculations are performed using the projec¬ 
tor augmented wave method implemented in the Vienna 
ab mitio simulation package (VASP) [50]. Perdew-Burke- 
Ernzerhof parametrization of the generalized gradient ap¬ 
proximation (GGA-PBE) is used for the exchange corre¬ 
lation potential [51, 52]. The plane wave energy cutoff 
is set to 300 eV, and the Brillouin zone is sampled by 
a 9 X 9 X 6 mesh. Based on the experimental or the 
optimized lattice structures, we apply the more sophisti¬ 
cated and accurate Heyd-Scuseria-Ernzerhof (HSE) hy¬ 
brid functional method [53] to the calculations of elec¬ 
tronic band structures, using a 300 eV wave energy cutoff 
and a 6 X 6 X 4 Brillouin zone mesh. 

In studying the surface states, we first use the bulk 
DFT results and the Wannier90 code [54-56] to construct 
the maximally localized Wannier functions for the p or¬ 
bitals of Bi and halogens. Based on the Wannier func¬ 
tions, we then use an iterative method [57] to obtain the 
surface Green’s functions of semi-infinite systems, from 
which we calculate the dispersions of surface states. 

Under strain, the atomic positions are allowed to 
be optimized until the force on every ion is less than 
0.01 eV/A by employing the vdW corrections [58, 59]. 
Since /3-Bi4Br4 has yet to be synthesized in experiment, 
we calculate the phonon spectrum to investigate its lat¬ 
tice stability, using the PHONOPY code [60]. The vdW 
corrections and lattice relaxations are also employed 
to obtain the optimized atomic positions and lattice 
parameters for /3-Bi4Br4, which are list in Table SI [21]. 
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